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Abstract. We prove that if t h-> u{t) e BV(]R) is the entropy solution to a N x N strictly hyperbolic 
system of conservation laws with genuinely nonlinear characteristic fields 

ut + f{u)x = 0, 

^— H then up to a countable set of times {tnlnsN the function u{t) is in SBV, i.e. its distributional derivative 

Uj; is a measure with no Cantorian part. 

The proof is based on the decomposition of Ux (t) into waves belonging to the characteristic families 

JV 

Ht) = J2 (*)^» (*) ' (*) e (IK) . fi (*) e ' 

r—K and the balance of the continuous/jump part of the measures Vi in regions bounded by characteristics. 

To this aim, a new interaction measure /lijump is introduced, controlling the creation of atoms in the 
measure Vi{t). 

The main argument of the proof is that for all t where the Cantorian part of Vi is not 0, either the 
I I Glimm functional has a downward jump, or there is a cancellation of waves or ttiG mGa,su.rc /^ijump 

positive. 
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1. Introduction 

In this paper we consider the entropy solution to the hyperbolic system in one space dimension 

iut + f{u)x = 0, u:R+ xR^ncR"", f e C^n,R^) 

\u{t = 0)=u ueBy{R;n), |t2|Bv<l- 

We assume that each characteristic field is either genuinely nonlinear or linearly degenerate, and in what 
follows we will refer to the unique solution constructed by vanishing viscosity or wave-front tracking, 
see [ETilO^ . 
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While linearly degenerate families do not gain any regularity during the time evolution, the genuinely 
nonlinear families show a regularizing effect due to the non linearity of the eigenvalue. The most famous 
one is probably the decay of positive waves, which in the scalar case iV = 1 takes the form 

nt < I, (1.2) 

where k is the genuinely nonlinearity constant, 

k = inf /"(m) > 0. 

u 

For a strictly hyperbolic system of conservation laws, even if the i-th family satisfies the genuinely 
nonlinearity condition 

DXi{u)n{u) > fc > 0, 

where {Xi}^i and {ri}^i are the eigenvalues and eigenvectors of A{u) :— Df{u) with a suitable ori- 
entation, then it may happen that new positive waves are created at a later time. In this case, the 



estimate (1.2 1 takes the form 



vntm<C^P^ + Q{s)-Q{t)}, 0<s<t, (1.3) 



s 

where i? is a Borel set, Q is the Glimm interaction potential and the constant C depends on k. The wave 



measures Vi{t) are defined precisely in Section 2.1 and roughly speaking they are the part of Ux{t) which 
has direction close to ri{u{t)) and travels with a speed close to Xi{u{t)). 

Since Q{t) is a decreasing function, an elementary argument yields that is absolutely continuous up 
to countably many times: in fact, if B is £^-negligible and v^{t){B) > 0, then by letting s t we obtain 

Q(t)-limQ(s)>^^^ffil>0, 

SO that t Q{t) has a jump downward. Being Q decreasing, this can happen only countably many times. 

A complementary estimate is the fact that also v~{t) has no Cantorian part. The first positive 
result has been given in [ADLj . where it is shown that the solution u{t) of a genuinely nonlinear scalar 
conservation law in one space dimension is SBV up to countably many times. In that paper, the authors 
consider the characteristic lines 

X = f'{u{t,x)), u{0,x)^y, 

and prove the following: every time a Cantorian part in Ux{t) appears, then there is a set of positive 
measure A such that all the characteristics starting from y € A are defined in the interval [0, t] but cannot 
be prolonged more than t. By the cr-finitncss of one can apply the same observation used to prove 
that the positive part of u^it) is abolutely continuous up to countably many times, and deduce that up 
to countably many times the solution u(t) is SBV. 

The use of the measure of the set A{t) of initial points for characteristics which can be prolonged 
up to time t has been applied to obtain extension of the above result: in |Rog| the SBV estimate is 
used for scalar balance laws, later extended to Temple systems in |AnN| and in [BDRj to the case of 
Hamilton- Jacobi equation in several space dimension with uniformly convex Hamiltonian. In the context 
of the Riemann problem for genuinely non-linear systems, the thesis has moreover been proved in |Daf] . 

The case of genuinely nonlinear systems of conservation laws is more complicate by the fact that 
centered rarefaction waves are created at t > 0, and thus the characteristics are not unique in the future. 
Thus, in estimating the /!^-measure of the initial points, one has to take into account also that interaction 
points can generate centered rarefaction waves, so that the estimate should be something like 

I initial points of characteristics arriving at t but not prolongable | 

> I measure of the Cantorian part of Wi(t)| - | amount of interaction in [0,t]]. (1.4) 

In this paper we use a different approach. Let Xi{t, x) be the Rankine-Hugoniot speed if u has a jump 
in the point (i, x) of the i-th family or the ?-eigenvalue of A(u) in the remaining cases. We first prove 
that not only the characteristic waves Vi{t) satisfy a balance equation of the form 

dtVi + dx{XiVi) = ^i. 
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with fii a measure bounded by the interaction-cancenation of waves, but also its atomic part Wijump(^) 
satisfy 



with /Zijump bounded measure. This measure fii 



differs from the interaction-cancellation measure 



because it is not when an atomic part in Vi is created, and it describes the natural behavior of solutions 
to genuinely nonlinear conservation laws: it is easy to create a shock because of the nonlinearity, but you 
can remove it only by means of cancellation or strong interactions. 

The second step is to use the two above balance equations to study the balance of Vi, Wi.jump and 



'^i,cont 



'^ijump 



in regions bounded by characteristics. The key estimate we obtain is that 



dtm > -c 



C\B) 

T-t 



M^^^({[t,r]xM}) , 0<t<r, 



(1.5) 



where the measure ^i^^'' is the interaction-cancellation measure ^i^^ plus morally the measure X^iLi iMiJumpI 
This is the companion estimate of (1.31, and using the same argument of the positive part we conclude 



with the main result of this paper (Corollary 3.3): 
Theorem 1.1. Let u(t) he the entropy solution of the Cauchy problem 



ut + f{u) 
u{t = 0) = 



0, 



u : . 



n c 



/ e C2(17;MN) 



for a strictly hyperbolic system of conservation laws where each characteristic field is genuinely non-linear, 
with initial datum u small in BV(M; f2). Then u{t) € SBV(M; Q) out of at most countably many times. 

1.1. Structure of the paper. The paper is organized as follows. 

In Section [2] the main notation and assumptions are introduced: strict hyperbolicity and characteristic 
families, and the decomposition into wave measures. A few fundamental results concerning the hyper- 



bolic systems of conservation laws are recalled: Lax's solution to the Riemann problem (Theorem 2.1) 
and Bressan's existence and uniqueness of a Lipschitz semigroup of solutions for small BV initial data 



(Theorem 1 2. 2 [) 

In Section [3| we prove Theorem |1.1| (Corollary 3.3 below). This result is a corollary of the fact that if 
the i-th family is genuinely nonlinear, then the i-th component Vi{t) of u^it) is SBV up to countably many 
times (Corollary 3.2), and the latter is a consequence of the estimates (1.3) and (1.5) (Theorem 3.1). 



In this section the notion of interaction-cancellation is recalled, and the SBV estimates are derived 
assuming (1.5), whose proof is postponed to Section [5j 

Since the proof depends on uniform estimates for the wave- front tracking approximations, in Section [4] 
we recall the basic properties of these approximated solutions. A key fact is the possibility to collect the 
jumps of the wave- front solution into two families: one is converging to the jump part of u^, and the 
other to the continuous part of u^,. This is done by defining the maximal (eq, £i)-shocks (Definition 4.1): 
< eo < £1 are two treshold parameters, fixing the minimal size of the jump (eg) and the lower bound for 
the maximal size of the jump (ei). This definition has already been used in jBre| to study the structure 
of the semigroup solution u, which we recall in Theorem |4.2[ From this result we obtain that the 
wave measure u,", vf 



2, jump 



converge weakly (Corollaries 4.5 |4.6 4.7). 

In the last section, 
distributions 



for the wave-front tracking approximation and the products , jump 



)n l5j we 



Section 5 we prove the decay estimate (1.5). First of all, we prove that the 



^t^ijump ~^ '9x(Aj l^ijump) 



are uniformly bounded measures: we denote them respectively by fj,^ and /ifjump (Proposition 5.1). The 



latter measure is called the jump balance measure. Since the bounds do not depend on the aprooximation 
parameter j/, it is possible to pass to the limit and to obtain the balance equations for the wave measures 
Vi, Wijump- The consequences of this fact however are not directly related with the SBV regularity, so 
we will address them in a forthcoming paper. Next, we study the balances of the measures , f^jump in 



regions bounded by minimal characteristics (Lemma 5.6). We then use an argument completely similar 
to the one used for the decay of positive waves in [Bre]: if / = [a,b] is an interval and Wj^^.^j^j (/) is 
too negative, then either v~^^^^p is cancelled or the characteristics a{t), b(t) starting from a, h collapse 
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in a future time; in the last case the time t for which a{t) — b{t) is of the order of the length of the 
interval divided by the amount of negative wave, or by the interaction measure, cancellation measure and 
jump wave balance measure in the region spanned in time by characteristics from the interval. In this 
way we give a precise meaning to the inequality (1.4 1. We thus obtain the estimate (1.51 first for the 
approximated wave- front tracking solution and for finitely many intervals, and then passing to the limit 
we recover the same estimate for the semigroup solution and for Borel sets (Lemmas 5.7 5.9 1. Adding 
the already known decay estimate for positive waves, we obtain the desired result (Corollary 5.111. 



2. General preliminary notions 



Consider the Cauchy problem (1.1). The following assumptions are done: 

(1) Strict hyperbolicity: we set A{u) — Df{u) and we assume that the eigenvalues {Xi}^ of A satisfy 

Ai(w) < • • • < An(u), u ^ ft. 
We denote the unit right eigenvectors, and the left ones satisfying r.i ■ Ij — dij, respectively by 

ri(u),...,rN(u) ?i(u), . . . ,/n(u)- 

(2) Each ^-th characteristic field is either genuinely non-linear, i.e. 

I DAj(u)rj(u)| > fc > 0, uefl, (2.1 : zGN) 

or linearly-degenerate, i.e. 

DX^(u)n{u) = 0, uen. (2.2 : zLD) 

By the general theory on hyperbolic systems of conservation laws, one has then the following theorem. 

Theorem 2.1. Let u — u^X{x<o} + W^Xx>a- Then there exists a unique self-similar weak solution 
whose shocks satisfy the Lax compatibility condition: 

Ai(w^) > Xi{uj^ ,uj^) > Ai(w^) at each jump [uj^,uj^] having speed Xiicu^ ,u}^). 

Proof. We recall just that one can define C^-curves \l/i((T)(uo) and the scalars Xi{uo, ^'^((7)1*0) by 

\f{M>,{a))-f{uo)^K{uo,^,{<7)){M>,ia)-uo) a<0' °' °' 



where is a vector parallel to and satisfying DAjf, = 1 if (2.1 : zGNl holds, otherwise |fj| = 1 
if instead (2.2 : iLD) holds. The proof (see e.g. Th. 5.3, |Brej ) is then based on the inverse function 
theorem applied to the local homeomorphism A : C/(0) 3 cr n- u+ G U(u^) 

A{ai, . . . ,crN)(u") = ^'n(o-n) o • • • o ^'i(cri)(u"). 

It is applied in a domain where the Jacobian of the map is uniformly bounded away from 0. From the C^- 
regularity of the curves its Jacobian at cr = is [fi| . . . [tn] (w^) = DA{a = 0)(u^). Under (2.1 : ?GN), 
if = ^j((t)(u~) with cr < we have a compressive t-shock, otherwise a centered t-raref action wave 
(which is continuous). In the case of linear degeneracy (2.2 : iLD ) we have an i-contact discontinuity. □ 

Theorem 2.2. There exists a closed domain T> C L^(M;r2) and a unique distributional solution u = 
u{t,x) = [u(t)](x) which is a Lipschitz semigroup T) x [0, -l-oo) — )■ T) and which for piecewise constant 
initial data coincides, for a small time, with the solution of the Cauchy problem obtained piecing together 
the standard entropy solutions of the Riemann problems. Moreover, it lives in the space o/BV functions. 

Below, we will refer to u as the semigroup solution, or equivalently as vanishing viscosity solution. 
When referring to pointwise values of u, we tacitly take its L^-representative such that the restriction 
map 1 1— > u{t) G L^{dx) is continuous from the right and u{t) is pointwise continuous from the right in x. 



SBV REGULARITY FOR GENUINELY NONLINEAR ID HYPERBOLIC SYSTEMS OF CONSERVATION LAWS 5 



2.1. Decomposition of into wave measures. By the smallness assumption on the BV norm of 
the initial datum, one can assume that the eigenvalues of A = D/ satisfy on 

inf Ai < sup Ai < ... < inf An < sup An, 

and the eigenvectors lie in different cones. We then decompose as in the literature (see e.g. pp. 93,201 
of [Bre| ) the vector-valued measures Ux, f{u)x along the right eigenvectors of ^ = D/. We adopt the 
following notation. 



(1) Under (2.1 : zGN) we normalize the z-th right eigenvector of A{u) so that DAjr^ — 1, as in the 



proof of Theorem 2.1 denote it by fi{u). Under (2.2 : iLDi just take fj(u) = so that 

(2) We fix the left eigenvectors . . . ,?n(u) so that li ■ fj — Sij still holds for i,j e {1, . . . ,N}. 

(3) Given two values G n,hy the solution to the Riemann problem, briefly recalled in Theorem 



2.1 



there exists a — (cti, . . . , (Tn) G such that u+ — A{(t){u ). We introduce the values 

u}Q^u-, = «'(cr,)(w,_i) i = 1,...,N. (2.3) 

We define Ai (m+ ,u^) as the i-th eigenvalue of the averaged matrix 

A(w,,Wi_i)= f A{eLj, + {i-0)Lj,_i)dd 

Jo 

and li (u"*" , M~ ) , Pi (m+ , u~ ) vectors satisfying li ■ fj — 6ij which are 



if (Tj < and (2.1 : zGN) holds, left\right i-eigenvectors of A^uj^^uj^^i) normalized so that 
Ai(wi) - A,(i:Ji_i) = k ■ (wi - a;j_i). (2.4a) 
• otherwise, li{u'^ ,u~) — li{0uJi + (1 — 0)LUi^i). Therj^ 

(7i = • (uj^ - oji-i). (2.4b) 
We call strength of the i-th wave the value 

C7i — k ■ {uJi — Wi-l). 

In the genuinely non-linear case, by the parameterization choice it is equal to Xi{ujf) — Aj:(a;j_i). 



Remark 2.3. Notice that, as in the proof of Theorem 2.1 we defined Ai(w+,w^) as the Rankine- 
Hugoniot speed of the z-th wave of the Riemann problem in case it is a shock or a 

contact discontinuity: indeed since 

f{uji) - /(wi_i) = / (D/(iLJi_i + 9{uji - uji^i)){uji - uji^i))d9 = A{uji,uji^i){uji - uJi-i), 







then, if the Rankine-Hugoniot condition f{uji)'-f(uji^i) — Xi(uji~uji^i){uji—uji-i) holds, Ui—Ui-i 
is an z-eigenvector of A{uji,uii-i) with eigenvalue \i(u^ ,u~). In case of rarefaction waves it is 
just an average speed, and we see below that in the wave front tracking construction one may 
choose this as the speed of the 'artificial' jump, up to small perturbations. 

We recall now the wave decomposition of Ux into wave measures, and consequently of the flow f{u)x- Let 

\i{t,x) = Xi(u{t,x^),u{t,x^)). 
k{t,x) = k{u{t,x^),u{t,x^)), fi(t,x) ^ ft[u{t,x^),u{t,x^)). 



"'^Denote by a{u^ ■,u~) the inverse function of cr i— > u+ = A{cr;u~). By tiie inverse function theorem J(7{u^;u~) = 
{DA) ~ ^ {cr{u~^ ,u~);u~). If either the i-th component ai is nonnegative or the i-th characteristic field is Unearly degenerate, 
by the ODE satisfied by 'I'i one finds Vcri(«+;'U~) = at least when m+ = ^i{ai){u~'), being in that case 'i/i{ai) = 

fii'^iicTi)). Then by the fundamental theorem of calculus and the definition of loq, - ■ ■ ,(^N 



Jo do 
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Lemma 2.4. Define the scalar measures Vi :— U ■ u^. Then the following decomposition holds 



N 

='^Vifi{t,x), (2.5a) 

i=l 
N 

ut = f{u)x = ^ Xi{t, x)v,fi{t, x). (2.5b) 

2 = 1 



Proof. The first equation is a direct consequence of the fact that fi , . . . , tn are hnearly independent and 
k ■ fj — 6ij. For justifying the second one, denote by the structure of BV functions (Sect. 3.7 of [AFPj l 

^ [Blu + Blu] +Y,Ht,lh{t)+) - uit,jh{t)-)]C\dt), 

where D^m,D!^u are the absolutely continuous and Cantor part of u^, while the sum is the jump part. 
The values u{t,jh{t)^) are indeed the values of the approximate jump at {t,jh{t)) by a fine property of 
semigroup solutions recalled e.g. in Theorem |4.2| below. Then by Volpert chain rule (Th. 3.99 of [AFP| ) 



/(«), = D/(u)[D^U + D^,li] +^ [/(u(i,7;,(t) + )) -/(«(i,7,(t)-))]£l(di).Dom(7.) 

N 

'^^(w)[D:u + D^,u] +^^[/(c.,(i,7.(t)+))-./(c..-l(t,7/.(0"))]^'(rfOLDom(7.) 

/iSN i=l 

N N N 

^\{t,x){v,)contri{t,x) + ^Xi{t,x){v,)j^^pr,{t,x) ^^\i{t,x)v,fi{t,x), 

i—1 i—1 i—1 

where we also applied in the last step that each fi{u) is a right eigenvector of A(u), and the Rankine- 
Hugoniot conditions of Remark |2.3| for the jump part. □ 

3. Main SBV regularity argument 



Given a semigroup solution u, the wave decomposition ( 2.5 1 of m^; along fi , . . . , f n reduces the vectorial 
problem to scalar ones. In order to check that there is no Cantor part in Ux, by (2.5) we will indeed check 
that there is no Cantor part in each Vi. 

The argument here is based on piecewise constant front-tracking approximations of the semigroup 
solution u, some recall is provided in Section |4j The magnitude of waves of each i^- front-tracking approx- 
imation may change in time only when discontinuity lines meet, and waves interact. In |Brel . Section 
7.6, an approximate conservation principle for wave strengths is expressed by introducing two finite mea- 
sures concentrated on interaction points of physical waves: the interaction and interaction- cancellation 
measures 

l^A{P}) = W^\ /^i. ({^}) = k cr I + S „ • (3-1) 

10 I ^ I 

Above a' , a" denote the incoming strengths of the two physical waves interacting at P. 
More precisely, given a polygonal region F with edges transversal to the waves it encounters (Figure [I]) , 
one considers the total amount Wj'j^*, W^out of positive and negative i-waves entering and exiting F. 
Considering the 'fluxes' (the amount of positive\negative ith-waves entering\ exiting F) 

in 

then one has the upper bounds 



W- = W- - W-^ W- , = - W'^, i = 1 N 

"m m "in out "^^out ^"^outJ ' -^7---7^'i 



m,,, - Wi^i < c/(F) \w^tt - wt\ < c/^(r). 

Since the above measures satisfy uniform estimates w.r.t. by compactness one can define measures 
/i^, /i^*^ on the entropy solution by taking w* -limits of the corresponding measures on a suitable sequence 
of i^-front-tracking approximations. 

Assume that the i-th characteristic field is genuinely non-linear, (2.1 : iGN). A balance of the kind 
above will lead in Lemma 5.9 below to the following estimate on a semigroup solution u: the density 
w.r.t. time of the continuous part (fi)cont of Vi — k ■ is controlled at time s by the Lebesgue measure 
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and by a finite measure on a horizontal strip around s. This estimate is analogous to the decay of positive 
waves proved in Theorem 10.3 of |Bre| . We state it here, postponing the proof to Section [5] 



Theorem 3.1. Suppose (2.1 : zGNl holds for the i-th characteristic field. Denoting 

{vi{t))cont + (%(0)jump : (%(0)jump purcly atomic, 
M+ X M such that for s > t > Q 



[v-,{t)\dt, v,{t) 
then there exists a finite, nonnegative Radon measure 



O 



{^)[c\B)/T + ^ll,^J{[s-T,S + T]xR)} > \{v,{s)Unt\{ 



B) 



V-B Borel subset of I 



(3.2 : i) 



The Radon measure fJ-i is defined in the next section as a w*-limit of measures (5.4) relative to 
front-tracking approximations. It takes into account both general interactions and cancellations and also 
balances for exclusively the jump part of v^. In the statement Ui(s) is the representative continuous from 
the right, coinciding with k ■ (u{s))x- 

Following an argument already in [ADLj . the estimate (3.2 : i) allows to prove that is made by 
merely a jump part concentrated on the graphs of countably many time-like Lipschitz curves and a 
part absolutely continuous w.r.t. the Lebesgue measure. Indeed, (%)cont is the integral w.r.t. C^{dt) of 
{vz{t))cont- if \{vi{s))cont\{B) > at some time s for an £^-negligible set B, then by (3.2 : «| the time 
marginal of fJ,^^'^ has an atom at s, which may happen only countably many times. 

Corollary 3.2. Let u be a semigroup solution of the Cauchy problem for the strictly hyperbolic sys- 
tem (1.1). Consider the t-wave measure = k ■ u^. If (2.1 : iGNl holds, then has no Cantor part. 



If all the characteristic fields are genuinely non-linear, by the wave decomposition (2.5) the above 
estimate yields then SBV([0, T] x M; M^) regularity of u for aU T > 0. 

Corollary 3.3. Let u be the semigroup solution of the Cauchy problem 

\ut + f{u)x = 0, 



u{t ^0)^u 



u : 



O C 



fec\n;i 



for a strictly hyperbolic system of conservation laws where each characteristic field is genuinely non-linear, 
with initial datum u small in BV(M; il). Then u{t) G SBV(M; fJ) out of at most countably many times. 



Proof of Corollaries \3.^ and \3.i^ By the general theory on 1-dimensional systems of conservation laws, 
u belongs to BV([0,r] x M; il) for all T > 0. By the structure theorem of BV functions (see Sect. 3.7 
of [AFPj ) . the derivative of u can be decomposed into a jump part, concentrated on a 1-countably 
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rectifiable set, a part which is absolutely continuous w.r.t. the Lebesgue measure, and a remaining part — 
the Cantor part — which is singular w.r.t. and vanishes on sets having finite HausdorfF 1-dimensional 
measure: 

Du Vu + D^^ It + DJ u. 

We want to prove that the Cantor parts u, D'^u of the both components Ut, of Du vanish. 

Denote by u{t) the 1-space-dimensional restriction of u at time t, which is a function of x. By the 
slicing theory of BV functions (Theorems 3.107-108 in |AFP| ). not only 

u{t)^dt, (3.3) 

but moreover the Cantor part of is the integral, w.r.t. C^{dt), of the Cantor part of u{t): one has also 
the disintegration 

D> = J {uit))dt. 

Since, moreover, by combining Volpert chain rule (Th. 3.99 of |AFP] ) and the conservation law one has 

D> = /'(u)D^,u, 

it suffices to show that for £-'^-a.e. t the function u{t) belongs to SBV(M;M^). 



If one sets Vi{t) = li-u{t)x, then by the slicing (3.3) and as in the wave decomposition (2.5) one derives 



Vi{t)dt i = l,...,N, u{t).x =^^f^{t,x)vi{t). 



1=1 



As a consequence, there is a Cantor part at time s in some Wj(s) precisely when u(s)j: has a Cantor part. 



It remains then to prove that (2.1 : zGN) implies that Vi{s) has no Cantor part for £^-a.e. s. By 



Theorem 3.1 the assumption of genuine nonlinearity (2.1 : IGN) implies that the estimate (3.2 : i\ holds 
for the i-th characteristic field. The fact that Wj(s) has a Cantor part means that there exists an C^- 
negligible compact set K with Vi(s){K) > and with no atom of Vt{s). Then for all s > t > 

0<\v,{s)\{K) = \{v,is)Unt\{K) 

0{1){£\K)/t + fi',^'\[s - r, s + r] X M)}. 
As we are taking C^(K) = 0, by outer regularity of Borel measures when t 4- this means that 

mP^(Mxm) >0. 

It is indeed true that if C^iK) = then ^J.(^^^i{s} x K) > 0{l)\vt{s)cont\iK) > 0. As the above measures 
are locally finite, this can thus happen at most countably many times: for all other times t the continuous 
part of Vi(t) is absolutely continuous. 



As a consequence, if ( |2.1 : IGN^ holds for all z" = 1, . . . , N then u{t) e SBV(M;M^) out of countably 
many times, the times when the time marginal of anyone of the various nj^'' has a jump, i = 1, . . . , N. 
This yields the membership of u in SBVioc(K+ x M; M^). □ 

4. Recalls on the approximation by front-tracking solutions 

We recall in this section a result about the convergence of a suitable sequence of j/-approximate front- 
tracking solutions (Pages 219-220 in [Bre| ) . 

The z/-approximate front-tracking solutions {u'^}i' are, roughly speaking, piecewise-constant functions 
obtained approximating by a step function the initial data, (approximatively) solving Riemann problems 
at discontinuity points and piecing together these solutions until the time they interact and discontinuity 
lines cross each other: at that time the procedure starts again. By the construction, which allows 
small perturbations of the speed, only two discontinuity lines are allowed to cross at one time. Each 
outgoing j-rarefaction wave [u~ , u"*" — ^',;(cr)(M^)] in the approximate solution of the Riemann problems 
is decomposed into small jumps [w/^w^^i = 'i'i{ah){uJh)] of strength ct/j at most i/; among the various 
possibilities, we let the h-th jump [ujh,u}h+i] travel with the mean speed Xi(u'^ ,u~). In order to control 
the number of discontinuity lines, if the interacting wave fronts are small enough a simplified Riemann 
solver is used, which leaves unchanged the size of the incoming waves introducing a non-physical wave 



SBV REGULARITY FOR GENUINELY NONLINEAR ID HYPERBOLIC SYSTEMS OF CONSERVATION LAWS 9 



front traveling with fixed speed higher then An; the total size of non-physical waves is controlled at each 
time by a constant ev 

At each time t, the restriction u'^{t) is a step function: its derivative consists of finitely many deltas. 
Below one can see that if one fixes suitable thresholds, it is possible to group these deltas in two families: 
up to subsequences, the largest of them converge in the i^-limit to the jumps of the entropy solution u 
at time for £^-a.e. time, while the others tail off up to the remaining continuous part of u[t)x. We 
distinguish the jumps, excluding interaction times, depending on their characteristic family. 
Looking at the (t, x)-variables instead of time-restrictions of u'', the derivative of u'^ is concentrated on 
polygonal lines and consists only of the jump part. Nevertheless, these broken lines can be grouped as 
follows in order to distinguish those converging to the jump set of u and those part of u'^ converging to 
the continuous part of Ux- 

Definition 4.1 (Maximal (eq, ei)-shock front). A maximal [eq, £i)-shock front of the i-th family of a 
v- approximate front-tracking solution u'^ is any maximal (w.r.t. inclusion) polygonal line (t,7(t)) in the 
(<, a;)-plane, with f < t < t^, 

- whose segments are i-shocks of u'^ with strength \(j\ > Sq, and at least once jfjl > £i > Eq; 

- whose nodes are interaction points of u'^; 

- which is on the left of any other polygonal line it intersects and having the above properties. 

The family of i-maximal (eq, £i)-shock fronts of u'^ is totally ordered by the displacements of the 
polygonal lines on K x K"*", and, up to extracting a subsequence, we are allowed to assume that its 
cardinality is a constant MJ ^ independent of v. Let 7!^'* s : ft^'''" n \ 1 — > M denote the 

uniformly Lipschitz paths of i-maximal (soj £i)-shock fronts in Ui, and consider their graphs 



™(eo,ei) 

'^(eo,ei)= U Graph (7('';^_^^) „^). 

m— 1 



Notice that the set J^f^ enlarges as (eoj^i) goes to 0. 

Consider sequences < < i 0, where £§ is not necessarily a power but has just an apex k. Up 
to subsequences and a diagonal argument, by a suitable labeling of the curves one can assume that for 
each i, /c, m fixed the Lipschitz curves 7!^'^ t,, converge uniformly 

^r4-?),™^^(-s,ef).™ asi/^00 (4.fa) 

to distinct Lipschitz curves which cover the jump set of u: out of countably many points in the {x^t)- 
plane, either u is continuous, and equal to the pointwise limit of u'^ at that point, or has jump. In the case 
of a jump, the jump point belongs to the graph of some 7* ^ for a suitable triple i, k,m; moreover, 

the left\right limits of the front-tracking approximation at the jump curve 7^^^ ^ converge, for m,k 
fixed and 1/ — ^ 00, to the left\right limits of u at that jump. Below one finds the precise statement. 

Theorem 4.2 (Th. 10.4 in |Bre) ) . The jump part of u is concentrated on the graphs of j'^^^k ^k-^ ^, 

N 

J=UUGraph(7(^,.^,.)J. 

i—l m,k 

Moreover, u is continuous and equal to the pointwise limit of u'^ out of J . Define the countable set 

e =60 U 61 U 62 U 63 

= I jump points of the initial datum ^1 1^ | atoms o//i^'^| 

IJ [intersection of any l\^k^^k-^ „^ and 7(jfe^gfc) ,„ with « 7^ j| (J [endpoints o/7(gfc ,,1-) „j|. 
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T/ie?|^ at each point (i,7(t)) = linii,(iiy, 7''(ty)) ofJ^\Q one has a jump between the distinct values 

u{t,-f{t)-)^limu''{t,,Yit.y)^ lim u(s,y), (4.1b) 

(s,y)^(t,7(t)) 
i/<7(s) 

u{t,-f{t)+) = \imu''{t^,j''{t^)+) = lim u{s,y). (4.1c) 

(s,y)-+(t,7(*)) 
■y>7(s) 

We now stress some consequences of the above analysis. On one hand, the fine approximation of u yields 
in turn piecewise-constant approximations of any function depending on time-space through sufficiently 
smooth composition with u, and depending on at jump points. The countably many points of O are 
not considered in these approximations statement, and not relevant. In particular, Theorem |4.2| yields 
the following corollary relative to the functions li, fi, A,; introduced at Pagejsj 

Notation 4.3. When adding an apex ly we refer to the corresponding quantities relative to the j/-front 
tracking approximation u'^: at each point Xi(t,x) := Xi{u{t, x'^),u{t, x")) and so on, while X'^{t,x) := 
Xi{u'^{t,x^),u'^(t,x~)) and similarly 

l'^{t,x) ■.= k{u''{t,x+),u''{t,x-)), r'^{t,x) ■.^n{u-'{t,x+),u''{t,x-)). (4.2a) 

Remark 4.4. Front-tracking approximations do not satisfy in general the precise Rankine-Hugoniot con- 
ditions, but they may have a small error. Then we define the functions, for i = 1, . . . , N, 

- \ Xi(u'^{t,x)) at points which are not on i-fronts of u"^ 

Aj (t, x) = < . . . . . . . . (4-2b) 

I 7*^ (t) at discontinuity points of the i-th characteristic family 

which are the actual speed of propagation of the outgoing jump of u'^ if belonging to the corresponding 
characteristic family. In general X'^ (t,^'^ (t)) ^ j'^it) and the difference vanishes as i/ — > oo. 

Corollary 4.5. At each point out of Q U the functions li, fi, Xi are the pointwise limits of their 
piecewise constant front-tracking approximations AJ', ZJ', fj' introduced above. 
Moreover, at each point {t,"/{t)) = liTa^{t^,j'^{ti,)) of J^\Q one has 

h{t,j(t)) = lim/r(t.,7"(t.)), n{t,j{t)) = limfni.,7"(t..)), 

A,(t,7(i)) = limAnt.,7''(i.)) = fini7'^(t). 

u u 

Proof. Out of O U J' the statement follows just by the pointwise convergence of u'^ to u provided by 



Theorem 4.2 because being / S C^(S1;K^) the eigenvalues\eigenvectors of A{u) depend continuously on 
the variable u. Also U, fi, Xi depend continuously on the values , being suitable means of continuous 
functions on the segment [uT , u^], and they depend on the time-space variables through the composition 
with the left and right limits u\u'^ at the jump: the claim at a jump point (t, ^{t)) = limi^(ty , ^^(tu)) out 
of 8 follows then by the convergence of u'^ {ty,^'^ [t^)^) to u^{t,"f{t)). □ 

Corollary 4.6. Each i-maximal {£q,£\)- shock front in Ui, converges (locally) as f oo to the cor- 
responding graph of j'^ ^ fc^ in the Hausdorff distance: (iij ( Graph7,^'fc , Graph 7* ^ fc, ) 0. 

Moreover, there is a pointwise convergence also for the slopes j'^^l ^{t) out of a countable set oft. 



Proof. This follows directly by ( |4.1[ ) and Corollary 4.5 taking into account that the v-bont tracking 
approximate solutions (approximativcly) satisfy on z-shocks the Rankine-Hugoniot conditions 

f{u+) - f{u-) ^ \{u+,u-){u+ - U-). □ 

On the other hand. Theorem |4.2| allows handy approximations of the jump part of Ux, and therefore 
of the jump part of the wave measures Vi. More precisely, consider for fixed (eo,£i) the measures 

"'jump --^(eoEi)' JU™P jump ' 



"^For obtaining the central term below we refer more precisely to formulas (10.78-79) in IBrel for the genuinely non-linear 
case, and to an adaptation of the proof after the absurd (10.83) in IBrel for the linearly-degenerate case. These formula 
show also that the curves in ( |4.1a| may intersect only at endpoints for k fixed, as we claimed before the theorem. For the 
purpose of this paper only the genuinely non-linear statement will be relevant, thus we have not specified that the curves 
7" are actually defined in a different way for linearly degenerate characteristics. 
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which are concentrated on finitely many segments, up to any finite time. 
Corollary 4.7. There exists a sequence Vk such that 



N 

^jump 

i=l 



DJ^ u = weak-lim ^ «^[ump ' (^'i)jump = weak-lim v-, 
Denoting by D^'^x — ^''^'^lR+xR absolutely continuous part, D!^w the Cantor part ofu^, one has 



N 

D^« + D^,« = weak-ilim«;:-l^-^^\ w^;!:^"^''^^ := (^,-^)^ _ ^ ...''-^■(^S,^?) 



jump 



Similarly each Xi{vi)jump is the limit of X^''v'^^^p^°'''^\ while \i{vi)cont the one of - t^j'um'p*'^'"^'''] • 

Remark 4.8. By inspection of the proof below, one can see that the convergence in the statement of 
Corollary 4.7 holds in the same way for the jump and continuous part of (^(i))^, for each time t except 
at most the ones corresponding to the points in 8 of Theorem |4.2| 

Proof. The BV structure of the semigroup solution u (Sect. 3.7 of [AFP] ) gives the formula for the jump 
part of Du; for Af^ = M^^^^, it is the matrix- valued measure 

N Ml 

D-i u = - (K) ri'H,^ = lim (m+ - M~) (g) nH^ / n. 

fc^oo^-^/ > LGraph{7', , ) 

z^l m— 1 ^ (eO'^l''"^ 

Moreover, each addend of the sum is the narrow limit of the analogous one relative to on the corre- 
sponding (eq, e5^)-shock front: by (4.1b I, (4.1c), Corollary 4.6 and the dominated convergence theorem 

Mi Ml 

(u+ - u-) ® nn^ I s weak-lim V 

^ ^ ' LGraph 7' >. y ^ ^ ' l 



m— 1 



Graph 7% , ly ^ ^ ' LGraph K"'!. , 



weak-lim -w"^-) ^n^-H^^^. . 



Collecting the two limits, one has the expression of the x-component 

N 



V)\ u = lim weak-lim ^ w 



k \ 

In particular, one can find a sequence Vk such that 



jump 
= 1 



D{. u — weak-lim 



N 

i/fc,i,(eS,ei) 



W 



jump 



i=l 



Since u'^'' converges to u in Lj'Q^(M+ x M), then the x-derivative converges weakly*. One then takes 
the difference between the two terms (u'^'')^ and w'!!' ' , and takes their weak*-limit, in order to see 

\ / -i- jump ' ' 

that each tt'cont weakly* -converge to the continuous part of Ux^ which is the difference between Ux 
itself and its jump part. 

The convergence of «;„^^'p^°'"^) = ^ -^j^'"'"'^ to («,)jump and of X-^v^^;"'''^ to A,(t;,)ju„,p follow as 
above, taking into account also Corollary |4. 5 1 □ 

5. Main estimate 
As explained in Section [3j the SBV regularity of u amounts to the estimate 

3C>0 c{c\B)/T + fii^''{[s-T,s + T]xR)}>\{v,{s)U^t\{B) V s > r > 0, 

which shows that in the case of genuine non-linearity the non-atomic part (ui(s))cont of Vi{s) is controlled 
by the Lebesgue 1-dimensional measure and the interaction-cancellation-jump wave balance measure of 
a strip around the time s. 
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In proving it wc consider the wave balance and jump wave balance measures 

^J.i ^ dtVi + d.j:{XiVi) Mi.jump = 9t(ui)jump + 9j;(AiZ;i)jump, 1 = 1,..., TV. (5.1) 



In Lemmas 5.3 5.4 below we prove that they are indeed Radon measures: they are controlled by using the 
interaction and interaction-cancellation measures /i^, /i^*^, in the i/- front-tracking approximation together 
with other terms vanishing in the limit. The negative part of ^ii^jump however may not be absolutely 
continuous w.r.t. fi^''^: this is why the statement holds with the measures /if'^'', which dominates ^^"-^ -|- 
iMijumpI for i = 1, . . . ,N. 

We will first prove estimates of the wave measures on piecewise constant ^/-approximate front-tracking 
solutions {u'^}] we then obtain the claims by passing to the limit on a suitable subsequence. 

5.1. Wave and jump wave balance measures. We consider the distributions 

^t{v^)jnmp + d^{XiVi)j^rnp =■■ l^Linmp Myump'^ ^* ("^''^ J,^'' J + ( A,<L j... ). 

Notice first that by Corollary |4.7| the distributions in the left column are limits of the ones in the right 
column along sequences (eq, Ei), as fc — > oo. This is why it seems natural to us to define them, even 
though we will apply in the present paper only the 'discrete' ones and we are not going to take advantage 
of the ones on the l.h.s.. 

For the rest of the paper we will often omit the index k: limits will be tacitly taken on subsequences 
of the one of Corollary 4.7 For every i/, {eq, Ei), i the following holds. 

Proposition 5.1. The distributions /i^, /iijump o.i^'d tJ-'( , t^'ifump'^^ '^''^ finite Radon measures. 

We prove in this section the proposition above, direct consequence of the more specific statements in 
Lemmas 5.3 5.4 We call these measures the i-th wave\jump wave balance measures respectively of u and 
u'^ . We introduce them with the aim of managing finer balances for the variation of waves, distinguishing 
different families and the part of their variation only due to shocks. This is technically more difficult 
and it is postponed to a future work: in the next section we give only rough balances on strips, for 
i^-front-tracking approximations, that will be applied in order to derive (3.2 : j|. 

Before stating the lemmas, let us illustrate in the scalar case what the two measures reduce to. 

Example 5.2. In the scalar case N = 1, the speed X{u) is just f'{u), genuine nonlinearity reads like 
f"{u) > fc > 0, the parameterization choice is 



I = 

The wave measure v = lux is then 



f"{u) at continuity points, 

— — _ — - at jump points. 



V = f"{u){I^lu + D» + ^'^""^l '^l^"'^ D> = (/'H)^, (o) 

where the last equality holds by Volpert chain rule (Th. 3.99 of |AFP| V 

In the case of a smooth solution one can directly compute that the two measures vanish: 

fi = dtv + dxiXv) = dt{f"{u)u^) + dx{f{u)J"{u)ux) 

= dt{dj'{u)) + = d^idtfiu)) + d,(nu)r{u)ux) 

= d,[f"{u){ut+f{u)ux)] =0. 
One can see that the measure if defined on the solution, vanishes also in the BV case: 



-J Vt{.f'iu))^ + J ipx{lut) = - J Vx{f'iu))t + J ^x{lut) 
^x[-{f{u))^ + {lu,)] ^^0, 
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where we applied repeatedly Volpert chain rule (for example at (•)) and at (*) we applied the conservation 
law -ut = {f{u))x = Xux- 

The measure fJ-jump defined on the solution, in the case of BV regularity, can be similarly computed as 

(°) /" ^ , ^ \^/ + -^T^i M T 



ump 



where we denoted by {7fc}fcgN Lipschitz curves covering the jump set of Ux- We now make some heuris- 
tics. Knowing that /i is a Radon measure, one obtains that the strength A'^(i,7fc(t)) — A^(i,7j.(t)) = 

lu-itjkit)) -^''^^^^^ ^^'^ ju^P 7fe is a function of bounded variation on time intervals where 7^ is 
separated from the other curves. In that case, if one had some regularity of u^{t,'^k{t)) one could derive 

mump = E { / Ht,jk{t))[{nu+)ut - nu-)ui) - A(/"(,i+x - r{u+)u-))]{t,j,{t))dt 

+ W[f{u+) ./'(^-)] }|(tO,,,(*0)) - Mf'{u+) /'(U-)] 

We now show that fi'l^ is a measure concentrated on interaction points, and mainly on interactions 
between physical waves, where it is controlled by the interaction measure /i^. At the interaction point 
{t,x), remembering that only one interaction may take place at one time, 

<(0(M) ~ hm<(t ~ e)(M) = n'^ {{{t,x)}). 

Lemma 5.3. The distribution /i^ is a Radon measure satisfying \^^ \ < 0{l)^l + for a purely atomic 
measure jp"! < 0{\)ei, concentrated on interaction points involving non-physical waves. In the w* -limit 
in V then one has j^J < 0(1)/!^. 

Proof. We proceed by direct computation of the distribution, that will be a measure concentrated at 
interaction points. Then we estimate the value at each point. 

1) Computing /i^ . Fix an index v and let {im}m"^i be time-parametrized curves whose graphs are re- 
spectively the discontinuity segments of u"^ . Given any test function ip G C^((0, +00) x M) one has 



L„ ■r,^ (I 

^ ' " {ip{t,e,n{t))y,-{u+ -U-)dt 



^Jr- dt 
in—l '"T- 

E 

m— 1 



(^(r+,C(T+)) - ¥^(r,;,C(T-J))[, 



As noticed in (2.4 1, if u+ = 5'i(o'i)u~ by construction 1^ ■ (u+ — u") = for j 7^ i, and for j i it is the 
strength of the wave, that we denote again by cTi (with an abuse of notation in case of genuinely non-linear 
shocks). The non- vanishing terms in the summation above therefore are the ones of discontinuities of the 
i-th family and non-physical fronts. If {(tj,, Zfe)}/; is the collection of 'nodes' of i-th fronts of for t > 0, 
the computation above yields 

^Pk5{r^,:,^) + p", (5.2a) 
k 

where, denoting at each node by a[, a'l the i-th component of the incoming strengths, cTi the outgoing, 

Pk = -g[- g'I (5.2b) 
while p" is concentrated on the nodes of non-physical waves, with 

P''[{P})=1] -o-r, -a', 

a = — u~, a' = u'^ — u'~ being the outgoing\ incoming strengths of the non-physical wave at P, if^ the 
outgoing\ incoming values of k on the non-physical front. Notice that a' or a" in the above expression 
may also vanish, for example in the last cr' = when a new non-physical wave originates. 
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2) Estimates . In order to bound /ij' one considers the Glimm functional of |Brej 

Tr" = I]ka|+Co ^ Co » 1, 

where a,/3 index the discontinuities at time i, with corresponding strengths cTq, ct^, and the second 
summation ranges over couples (a, (3) whose discontinuity points Xa < xp either belong to families 
ka > ki3 € {1, . . . ,N + 1} or belong to a same family but at least one is a genuinely non-linear shock. 
Let T'' be the same functional, but neglecting in both the summations all terms involving strengths of 
non-physical waves. Since the total strength of non-physical waves at fixed time is controlled by Ei, by 
construction of ^/-approximate front-tracking solution, then — T'^\ < 0{l)e^. 

By choosing Co big enough and the smallness of the total variation, by those estimates in Lemma 7.2, 
Page 133, of |Bre| and Pages 137, 138, at each interaction at time r 

< 0{l)\a'a"\ < 0(l){Y^(r,7) - T''(r+)}. 

If %P is the set of interaction times involving only physical waves, by definition of the interaction measure 

|/^rk-XPxM< 0{l)^li^^.^^, fii{R+ X E) < T''(0). 

For the set T" of interaction times involving some non-physical front, instead one has the bound 

Kl< E 0(l){T'^(r-)-T''(r+)} 

J2 C'(l){T''(T-)-T'^(r+)-[T'^(r-)-T'^(T+)]} 

< Oil){r-^ir-)-r'^{T+)-[T-^{T-)-T'^{T+)]}+0{l)eM\ 

- O(l){[T''(0) - T^(0)] - [T^icx,) - T'^M]} + Oil)eM < 0{l)s,, 

because at interactions involving also non-physical waves T'^ does not vary (o), while at interactions of 
two physical waves a' , a" e.g. of the i-th family most of the terms cancel (★): 

X%T-) - T''(r+) - {T\t-) - T'^(r+)] < 0(1) ("^ |a'^| + |a' + a" - a\\ ^ \a,\. 

k^i (3 non physical 

One can gain the estimates on fii in the z^-limit because T''(0) can be bounded uniformly in ly. □ 

The following lemma deals instead with the jump-wave balance measure, which is again concentrated 
on interactions. At the interaction point (i, a;) one has 

<mieU)) - - ^i^^ie^d - ^^({(^' ^)})- 

The positive part is absolutely continuous w.r.t. the interaction-cancellation measure, up to a remainder. 
The bound on the negative part is obtained instead considering the history of the interacting waves, and 
involves the Glimm functional T. 

Notations are illustrated in Figure[2] Denote by {{ti^, x^i)}^ the finite set of terminal points of the i-th 
(eo,£i)-jump set ^(e'Jgj) of it" and let {'yh}h, {f/ij/i denote the relative maximal (eo, £i)-shock front and 
the strength < Eq just after becoming below the threshold (if any, otherwise vanishes). 

Lemma 5.4. The positive part of satisfies the following inequality 

mS^p"^ < Oil)f^i'^' -Y^auS^u^,,) \au\ < 0(1)-^/.^^ (7.)- 

h 

Also the negative part o/Mrjump^^ concentrated on the set of nodes of J^^f^^ and its mass is uniformly 
bounded by the Glimm functional T. In the w* -limit of Corollary \4. 7| one gets 

-P < Aiijump < 0{l)fj.^'^ p>0 finite. 

Proof. As before, we proceed by direct computation of the measure. 
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Figure 2. Discontinuity lines of u'^ belonging to the z-th (eq, ei)-jump set J^'^'J^^^. 
Dashed lines either are below the threshold (e.g. cti) or belong to a different family j ^ i. 



1) Computing A'^j^u^p^'' • One has to take into account that some of the shocks of u'^ are not present in 
the maximal {eq, ei)-shock fronts, and none of the rarefaction fronts and non-physical waves appears. Fix 
an index v and let {jm}m=i be time-parametrized curves whose graphs are respectively the discontinuity 
segments of u'^ present in >7(^^^ ^^y Given any test function ip G C'^^{{0,+oo) x M), and omitting for 
simplicity most of the indices v, (eq, £i), * from the second line below up to the end of the proof, one has 

~-^(eo,si) 



= XI (¥'(^m>7m(T-,i)) - 'P(T-„,7(T-,;)))cr^. 

m 

If {(Tfe, Zk)}k are the nodes in jT^'^'J for t > 0, then the computation above can be rewritten as 

Ihinl^P^ = 1k^(rk ,z,), (5.3a) 
fc 

where, denoting at each node by a', cr" the i-incoming strengths, a the outgoing, one has 



a — a' — cr" at triple points of JT^'^'J 

a at initial points of maximal fronts of J^^^^^ ^^y 

— cr' when a maximal front of JZ^^'J ends without merging into another, 

a — a' otherwise. 



(5.3b) 



2) Upper bound. At triple points of J^^^^ one has \qk\ < 0{l)fj,l, by the interaction estimates in 



Lemma 7.2 of [Brej and by definition of the interaction measure (3.1 ). 
At initial points of iT^^'* ^^y by genuine nonlinearity qk < 0. 
At internal nodal points of a front where 

- another shock of the same family, not belonging to any front in JT^"^'* ^_^y merges, again qk < 0; 

- a rarefaction wave front interacts, a cancellation occurs and \qk\ = |cr — cr'| < filf^ {tk,Xk) by 



definition of the interaction-cacellation measure in (3.1) 



any curve of different family interact, interaction takes place and \qk\ = |cr — cr'| < 0{l)^l{tk,Xk) 
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At terminal points, since the shock gets cancelled from Jj^f^ there must be a cancellation by a 
rarefaction front or by a wave of different family. If 7^ is the maximal (eq, £i)-shock front through the 
point, by inspection one can see that the strength must increase from a value (Tq < — £1 at some point 
up to a value < — Eq on the last segment, which colliding the other wave front becomes cti and exceeds 
— eo- By the interaction estimates, at nodal points of each front the strength may increase at most of the 
amount of interaction-cancellation at that point. One can then conclude 

£1 - eo < kol - |ai| < Tot.Var.+(u,7fe) < 0(l)/ii^(7fc). 

This yields the final estimate, at terminal points: 

qk = -<yi + {<yi - Qk) < -^^{ei - £0) + 0{i)fii{tk,xk) < -^^0{i)fii^{-/k) + 0{l)^li{tk,xk). 

£1 — £0 £1 — £0 

Notice that cri is one of the ak in the first statement. 
By summing up the different contributions, being 

£1 — £0 

one gains that /ij'j^u^p^^ ~ ^{^)^^^^ + J2 ^^hSit^^xh} is a signed measure, and thus ^J■if^mp^^ ^ Radon 
measure. We got an upper bound, uniform w.r.t v, of its positive mass, and we also got that the limit 
in V in the sense of distributions, i.e. /iyump, is a Radon measure. Notice that ^ |(T/,.| by the estimate 

'0- 



above vanishes in the limit of Corollary |4.7[ as we have chosen £^ > 2*''£q 



3) Lower bound. We derive now a lower bound uniform w.r.t v of the negative mass of ^J'ij^°J^p^\ suitable 
to estimate that the negative part of the limit /Xijump is actually a finite measure. 

Consider the nonnegative measure p. = — Mr/ump^^ + + J2 \^h\hth,Xh)- Consider for a > the 

Lipschitz test function (pait) = X[o.T+a]{'t) — (t — T) /axiT.T+ajit)- being the time marginal v'jl^^^p(t) of 
wj'jjjjjjp absolutely continuous and locally finite one has 

/2([0,T] XR)< f ^^dfl = - f ^a<jump'^ + I Vc.df^i'' + J2 



< 



h 

((^Ji + Ar((/.o).)d[<j,,„p(t)]dt + [<jump(0)] (K) + ([0, T + a] X 



T+a 
T 



{t)] {R)dt + Kj,^p(0)] (M) + O(l)Mi^([0, T + a] X 



1 

a 

By the uniform bound on ji'f ju,„p(^)|j as a J, one obtains 

//([0,r] X M) < Ao^^Tot.Var.(<j,,,p) + (!^(1V^^([0, T + a] x M) < O(l)T^(0). □ 

5.2. Balances on characteristic regions. We already discussed the distinction between the part of 
which approximates the jump part of Ux and what remains: now we denote shortly 

^i,cont ^2, jump Mi.cont Mi Mzjump" 

In the proof of Lemmas |5.3[ |5.4| we derived the balances 

<(i) - = ^^^m X M), <jump(0 - <jump(i") = <jump({0 X M). 

We estimate in this section balances of these measures in regions bounded by generalized I-th characteris- 
tics, based on the assumption of genuine non-linearity of the i-th characteristic field. A relevant measure 
will be the approximated interaction-cancellation-jump balance measure 

M-'^-'^M^'^ + lM-rjumpI- (5.4) 

The measure //j'^'' is defined as a w*-limit of the measures ^J^i'x,''' ■ 

We recall that a (generalized) i-th characteristic is an absolutely continuous curve x{t) satisfying for 
a.e. t the differential inclusion 

xit) e [X.,{u{t,x{t) + )),X,{u{t,x{t)-))]. 
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Due to the presence of discontinuities, if u is a semigroup solution of a system of conservation laws 
or a i/-approximate front-tracking solution there are several backward and forward characteristics, and 
they may collapse. They are polygonal lines whose direction changes at interaction points or hitting 
discontinuities, which are their nodes. 

One can select the minimal i-th characteristics of the piecewise constant approximations starting at 
to from any point x, which means 

x{t;to, xq) = mini x{t) : x{t = to) ^ xq, x{t) E [Xi{u{t,x+)),Xi{u{t,x^))] 



Analogously, one can select the maximal one. 

For the rest of the paper we will select for simplicity minimal generalized characteristics. We stress 
however that the same statements hold if one selects any other family y'^{t;to,x) of generalized charac- 
teristics of u'^ which, as the minimal and maximal ones, have the semigroup property 

y^ik; to, xo) ^xo; y''{t + t; to, xo) = y^ir; t, y^it, to,Xo))- 

Notation 5.5. Given an interval / = [a, 6], we define the region bounded by the minimal i- 

characteristics a{t),b(t) starting at to respectively from a, 6, and its time-section I{t), as 

A*,",'l]^{it,x): to<t<to + T, a{t)<x<bit)}, /(<) = [«(*), 6(t)]. (5.5) 



[a.b] 

M t- 



Given J — IiU ■ ■ ■ U I^i the union of any disjoint closed intervals /i, . . . , Im, M e N, we denote by 



the union of the regions A °^'' . Similarly, J{t) = A}""" n {i} x M = U • • • U /m(0- 



Lemma 5.6 (Approximate wave balances). Assume (2.1 : iGNl and consider regions hounded by i- 
characteristics as above. Then for any T,to > and t — to + tq 

<(i)(j(t)) -<(io)(J) < Oil){^,'f{A'r)+e.), 

<co„tW(^W) -<co„t(io)(j) < 0{l){fi',^/{A'r)+s,). 

The lemma shows that the variation of both the 'continuous' and 'jump' part of the wave measures on 
characteristic regions is controlled by finite, purely atomic measures, concentrated on interaction points. 
While for points in the interior of the region a perfect balance holds, with the suitable measures, in the 
proof one defines 'fluxes through the boundaries' in order to take into account the variation due to waves 
that enter\exit the region. Then the fluxes are controlled by the assumption of genuine non-linearity and 
the interaction-cancellation measure, by the choice of the speed of the rarefaction as a mean of the values 
A(u±). 

Proof. Fix for simplicity of notations to = 0, AZ_ = Aj^. Given a closed interval [a,b] we briefly denote 

<(i):=<(t)([a(t),6(t)]). (5.6a) 

We prove before the first balance, which is most of the work, and then we explain the others. 

0) Contribution of non-physical waves. The total strength of non-physical fronts at each finite time 

is controlled by e^, as well as the mass \p'^\ of /x"^ due to interactions involving non-physical waves 



(Lemma 5.3). Since thus non-physical interactions and waves are not relevant, they may be neglected 
by removing from both ^'^ , the relative terms: denote by p,'^ = ^i'^ — p'^, v'^ = w''L{physicai waves} this 
simplified measures. It suffices to obtain the claim for these measures, then the other follows. 
1) Instantaneous estimate. The region f,] is bounded by generalized i-characteristics: then the 
function v^{t) may vary only at times t where an interaction takes place in [a{t),b(t)], an i-rarefaction 
wave leaves the region AJ^ ^ or an i-shock enters. The following cases may occur: if physical waves in 

A'^^ change at (i,x), then one can define 'I'j'™^] by 

<(£+) - <(£-) = p^iliix)}) + ci>_:;-^,({(£,f)}). 

Formula ( |5.2[ ) and the inspection of all the possible cases, that we perform without taking into account 
that a{t),b{t) are minimal, yield that ({(i, i)}) is different from only if x coincides with a{i) 

or b{t). The value can be computed exactly. The most relevant cases are described below, and they 
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correspond to interaction among physical waves of the same family. The others are analogous. For 
example, we do not consider when only part of a rarefaction exits the region: one can see that the value 
of is controlled by above and below by two cases we consider (the whole rarefaction exists, or no 

part of it exits) . In these cases, that we illustrate in Figure 3 ^^'iTm is precisely equal to 





— a; cr > If an «-rarefaction wave exits. 

c'; a' < If an i-shock enters. 

cr' + cr"; cr', cr" < If two i-shocks out of A^^ ^ interact. 

cr'; cr' < If two i-shocks, one in one out of ^[a,fc]' interact. 

If an i-shock out of AJ^ j,j interacts with an i-rarefaction wave, 
and a «-rarefaction wave exits the region AJ^ . 

If an i-shock out of AJ^ j,] interacts with an z-rarefaction wave, 
and no i-rarefaction wave exits the region AJ^ ^ . 



-a + a'; cr' < 0, ct > 
cr'; cr' < 



If an i-rarefaction wave out of AJ^ interacts with a i-shock in 
AJ^ , and a i-rarefaction wave exits the region AJ^ . 

If an z-rarefaction wave out of AT„ ,,i interacts with a i-shock in 
cr > - • ' J 

— ^fa 6] ' ^^"i ^'^ i-rarefaction wave exits the region AJ^ . 

If an z-rarefaction wave and i-shock both out of AL interact, 

f a"; (t' > 0, ct" < , _ , ■, ^.x, ■ Ar 

' — ' — and no z-rarefaction wave exits the region A j^ ,,]. 

If an z-rarefaction wave and i-shock both out of ^[a,&] interact, 
and a i-rarefaction wave exits the region AJ^ . 



cr"; CT,cr' > 0, cr" < 



The explanation is that when a physical interaction takes place at the boundary one has to balance the 
part of the waves that are out of the region we are looking at. The assumption of genuine-nonlinearity of 
the i-th characteristic field is what determines the sign of the big discontinuities, the i-shocks. Considering 
that, the contribution of may be positive only in the last four cases, when an z-rarefaction wave 

comes from the outer region and hits a shock: in this case we now show that cancellation occurs and $j 
is controlled by the amount of cancellation in /i^*^ at that point. Indeed, suppose that an i-rarefaction 
£' hits an z-shock £": the value Um of between the two fronts should be the same, while the incidence 
condition implies i' > £" . Moreover the speeds i',£" can be estimated by the intermediate value Um and 
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the strength of the waves: one has 



-3(7" > a' 



This says that if a cancellation occurs, then the shock should be of the same size as the rarefaction. 
Therefore the rarefaction even when not cancelled is controlled by the amount of cancellation. Then 

More formally, the 'flux' $j could be defined at once considering the exterior trace of u'^ on AJ^ . 

2) Single interval estimate. By the choice of AJ^^, which contains /(r) but not /(O), there is no 
flux concerning waves exiting at r and no one concerning waves entering at — even if an interaction 
takes place at r or at — because Vi{t) is continuous from the right. Indeed interactions at t = are not 
taken into account by /ij on AJ^ ^^j. We recall that now we temporary neglect non-physical waves. 

Piecing together the various physical interactions one obtains 

v^r) - ^r(O) = f^U^laM) + ^u:k^GTaMa) U Graph(6)) ^'^'^^^^{B) < (5.6b) 

where B C Graph(a) U Graph(5). This formula holds also in the case a ^ b. When removing the hat 
the same equation holds with an additional term uniformly controlled by e^, however we prefer to take 
non-physical wave into account after having obtained the balance for countably many intervals. 

3) Estimate on more intervals. Consider now two disjoint closed intervals Iq — [xq,Xi], Ii = [2:2, CC3] 
and denote the relative selected i-characteristics starting out from by Xh{t), Ih{t) = [x2h{t),X2h+i{t)\- 



Notice that the same equation as in (5.6b I holds immediately for Aj^^jj^ on time intervals where a;i(io + 
r) < ^2(^0 + ■'")■ We consider thus the case xi{t) = X2(t) at time t < Iq + t and then, by the semigroup 
property, up to time tg + t. 



JO+T 




Figure 4. Illustration of Lemma 5.6 balance on two intervals 



As illustrated in Figure 4l decompose A*,"'^ U A*,"'^ into the regions At"'^ ,^1^°'^ ,^f'^/-T ^-m, with 

l_J ^ JO Jl ° Jo ' Jl ' [xo(t),X3(t)]^ 

t' = t ~ to, where one applies separately (5.6b). In order to shorten notations, fix again to = 0: 



v-^{T){[xo{r),x,{T)])-v!^{0){loUh) 

= v-^(r){[xo{r),xs(r)])~v-^{i){[xo{i),xM) 

+ v^{i){[xoii),x,(i)])+v;^ii){[x2(t),x,(i)])-v^ii){{^^^^^ 

-t)r(o)(/o)-t)r(o)(/i). 

If one denotes briefly $j ((i, t]) — ^^'^^"^ {J^il^ and so on, the balances yield 
t)r(r)([xo(r),X3(r)])-z;r(0)(/oU/i) 

- (^;r„^1.3(,-)j) + Ar(4) + Ar (4) - ^^m^^m 
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It remains then to notice that as the regions Aj^, overlap at the interaction point then 

Ar(<^,.3(f)]) + Ar(4j + Ar(4) = Ar(4u/J +Mr({(Mi(^))}). 

Moreover, the waves entering at (i,xi(i)) from outside the are counted both in 

$^'^"'((0, i\) and <i>j j"'((0, t]), but only once in ((*-*' '''])' while the ones interacting at (£, xi{i)) and 

coming from the interior of ^/gu/i ^^"'^ counted precisely once in 4>^'™'((0,t|) + <J>^'°^"*((0, since the 
amount of these waves from the region in between is precisely — p.!^{{{i,xi{i))}), therefore 

Collecting the terms one obatins 

<(r)([xo(r),X3(r)]) - «r(0)(/o U h) = fi^ {A^ U A^) + ^^'ouV ^D- 

The reasoning applied to two intervals Iq, Ii can be applied as well to any finite number, obtaining the 
analogous equation. Since the flux can be estimated by the interaction-cancellation measure, one finds 
then that for all M G N and real intervals /i , . . . , Im 

v^{T){h{T) U • • • U lM{r)) ~ <(0)(/i U • • • U /m) < (Ar + 0{l)f,'f){Al U • • • U A} J. 

Since the total strength of non-physical fronts at each finite time is controlled by e„, as well as the 
mass \p'^\ of /ij due to interactions involving non-physical waves and fj,^ (Lemma |5.3[), from above 



t;r(r)(/i(r) U • • • U /mM) - <(0)(/i U • • • U Im) < {^^^ + Ku-u/.J + 0{l)e. (5.7) 

<0{l)[t/^'{Al^...^,J+e,). 



4) Balance for the jump part. For the 'jump part' by (5.3) one can repeat the same argument as 
above with the relative measures, but the fluxes are non-positive because only j-shocks are involved, so 
that one does not need to estimate the positive part (which vanishes) of the fluxes by the interaction- 
cancellation measure: therefore 

<iun.pir){HT),b{r)]) - <j,^p(0)([a, b]) = Mr,jump(^[a.^) + *-::Mr"'(^[a,b]) (5-8) 

holds with new fluxes (i)':''°"*'J"™P which take into account only the contribute in ^'^f^t due to z-shocks 
in JT^'^'* In particular this flux is non-positive. Explicitly the new flux at the boundary is 

If an i-shock of strength a' in j7/'* n enters, possibly interacting with a front 
a'; a' <0 ^ ^""'^^ .1 

wave either in AT , , or not belonging io J, ' v 

a' -\- cr"; fj', a" < If two j-shocks out of AJ^^ ^ and in J^^^^ interact. 

Otherwise. 

The generalization to countably many intervals holds precisely as before. 

5) Balsince for the continuous part. By subtracting (5.8) to (5.6b) one finds the analogous estimate 
for Vcont, which again holds in the same way for countably many intervals, with a new flux at the boundary 
^i'[ah\°^ (^fa b]) '^hich is the difference — ^"'[ab]^'^'^ '■ ^^otice indeed that one can take advantage 
of'°' 

out _ out jump ^ r^,. out _ out jumpn+ ^ T^''^™* 1+ < „/C^ 

Indeed, the negative terms in '1?^ cancel positive ones in , without adding a positive contribu- 

tion in the difference. The last equation in the statement holds then by definition of nl^'^ — M^,*" + |A*ijump I ' 
which controls - Muump- ^ 
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5.3. The decay estimate. We prove in this section an estimate analogous to the decay of positive 
waves (Pages 210-216 in |Brej ). but for the negative part of ('Uj)cont, under the assumption that the z- 
characteristic field is genuinely non-linear. We bound from below the continuous part of the waves of a 
semigroup solution u on a Borel set B at time t by the wave\jump wave balance measures fJ^^\^i^^jump on 
a strip around t of arbitrary height t < t and the Lebesgue measure of B divided by r: together with 
the previous result due to Bressan we obtain 

IWcontlW(S) <(^(l)'C\B)/r + 0(l)(/4^'^)([t-r,i + T] xM) > r > 0. 

We remind that we defined /ij*^'^ as the i(;*-limit of the measures ^l'-^'^ in (5.4). 

We remark that here there is no sharpness purpose. We first prove an estimate from below for the 
'absolutely continuous' part of w^. As a consequence of Corollary 4.7 and Remark |4.8[ when passing to 



the limit we get the claim (Lemma 5.9). 

Lemma 5.7 (Approximate decay estimate on intervals). Assume the i-th characteristic field is genuinely 
non-linear. Then for any disjoint closed intervals {Ih}heN '^''^^ r, > one has the bound 

-<cont(io)(/i U • • • U /m) < 0{l)[c\h U • • • U lM)/r + li^,' {A]%,„^,J + + £i}, 

where ^/°u...u/jvj region hounded by generalized i-characteristics in (5.5). 

Remark 5.8. Under the assumptions of the above lemma, an analogous statement holds for any Borel 
set B if one considers its evolution by minimal generalized i-th characteristics (or a different selection 
with the semigroup property). Denoting the generalized z-th characteristics {y'^{t; to, x)}x£ b and defining 

Ab = {2/"'(i;io,a;)}?=^^<^<t„+^ one has 



e„ + eij. 

Indeed, this holds true for countably many intervals by Lemma |5.7| and the observation that 

i^'3'i^l::l]) = M(^^(^[a,bi) = ii',!i-'(Ai;j]n{t > to}), 

for any choice of characteristics starting out from [a, b] and having the semigroup property. Since on the 
r.h.s. we have a nonnegative Radon measure, by inner/outer regularity it can be extended to Borel sets. 

Proof of Lemma \5. 7| We adapt the argument in [Bre] (§10.2) about the decay of positive waves, whose 
main steps we briefiy recall here. Since the generalization to countably many intervals is analogous to 
the one of Lemma 5.6 we consider the case of a single interval / = [a, b]. 

The heuristic is the following. Suppose w|^cont(^o)(^) < and there are only i-fronts in /, which is the 
main case we would like to manage. Then initially C^{I{t)) decreases at a rate at least (to)(^)/4; if 
C^{I{t)) keeps on decreasing at least of that rate between times to, to + one can estimate the initial 
value C^{I) from below as stated just by integrating this differential relation. Otherwise, interactions 
must take place in order to decrease the rate, with a special care for the ones at the boundary. 

For simplicity of notation, we set to = and we consider the «-th wave measures omitting the indices 
I. Let the interval [a, 6] evolve by the minimal forward characteristics: its length at time t is 

z(t) = b{t) - a(t). 

This function is absolutely continuous and it satisfies z{t) = A(t, 6(t)) — A(t,a(t)) for a.e. t. 

Estimate on the influence of the other families on the i-wave measure. At Page 213 of |Bre| it is 
explicitly defined a function $ which contains the variation of speed due to the presence of waves of other 
characteristic families. $ is piecewise Lipschitz continuous, non-decreasing, with discontinuities only at 
interaction times. More precisely one has the estimate 

iz(t) + e(t) - v^{t)\ < o{i){e. + mm) (5.9) 

where v'^{t) counts the variation of speed due to the waves of the same family, and ^ corrects it taking 
into account waves at the boundary: we adopt the notation 

<co„t)W :=<co„t)W([aW,K^)]), 
m ■■= (A(t, ait)) - A(t, ait)-)) + (A(t, fe(t)+) - A(t, 6(t))). 
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Case 1. If z{t) - ^{t)z{t) < <Q„t(0)/4 for all t, then 

d/dt{e-^o'^z{t)) = e-^o^lzit) - 'i'{t)z{t)} < e^-^'o *w'^cont(0)/4 < w''cont(0)/4 

and integrating the inequality between times and t, since z{t) > 0, one has 

-^(0) < e-fo-^zir) - z{0) < r<,„,(0)/4. 

Case 2. Suppose instead z{t) — <i>(f)z(f) > wJfQjj^(0)/4 at some time t. 
By (5.9) we have that 

z{t) - ^{t)z{t) < v^it) - ^[t) + 0{l)e,. 

Observing that \(t,a{t)) is a suitable mean of \{t,a{t)^) and vy^^-^^{(t,a{t))) — A(t,a(t)+) — X{t,a{ty 
by (2.4a), one derives the bound 

> <umpW-2ei. 
The two inequality, being v^^^-^^ negative, yield immediately that 

m - mz{t) < + i,f„„,p(0 - m + 0{l)e, < <,„t(i) + 0{l)e, + 2ei. 

At time i, taking also into account that we are assuming z{t) — <P{t)z{t) > t>^Q,j^ (0) /4, we get thus 

<„„t(0)/4 < <o„t(i) + 0{l)e. + 2ei 



By the balance for the continuous part in the proof of Lemma 5.6 we obtain 
Conclusion. Collecting the two cases one finds the bound uniform in i> 



<„„,(0) > - a)/T + fii'^'iAl,^) + £1 + e.j. 



This shows the claim for a single closed interval. Similarly to what is done in Lemma 5.6 the argument 
extends by direct computation to the union of countably many closed intervals. □ 

We finally arrive to the decay estimate for the negative part of (wi)cont- The characteristics of the semi- 
group solution u we get in the limit in general are no more minimal, but just generalized ^-characteristics: 
this is why in (3.2 : i) we have chosen to state it enlarging the set As up to an horizontal strip. 

Lemma 5.9 (Decay estimate on intervals). Assume the i-th characteristic field is genuinely non-linear. 
Then there is a choice of generalized i-th characteristics {y(i; s, a;)}^gR of u such that for any Borel set 
i? C K and t^ > t > Q the following estimate of the i-th wave holds 

-{y,)cont{to){B) < 0{l){c\B)/T + f,','''(A^)} Ab = {y{t-to:x)r,fZ<t,+r- 

Remark 5.10. By the Jordan decomposition of (fj)cont(io) ^^id the arbitrariness of _B, this is precisely a 
lower bound for the negative part of (f j)cont(io)- 

Proof. As usual, simplify notations by setting t^ — and omitting i. 

When io was not a time in Q of Theorem 4.2 by Remark 4.8 the 'continuous part' of v^{Q) w*-converge 
to the continuous part of v . Nevertheless, we do not know whether the positive and negative part [w''(0)]* 
of v^{Q), and [i'f,cont(0)]'''' converge to the positive and negative part in the Jordan decomposition of w(0), 
and (w(0))cont- However, if one defines measures [v\^ by the relation 

[<o„t(0)]^ ^ [vf. 

where the convergence is obtained by compactness and up to a subsequence, then 

[V]+ ~[v]- ^{v{Q)Unt > [(^^)co„t(0)]+ > [(v)cont(O)]-. 

In general \v\^ are not orthogonal. 
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Since the atoms of [v]^ are at most countably many, one can consider first finitely many intervals 



whose boundary has \i 
/i U • • • U Im converge: 

]^(0)(J) 



-measure, where 



Then the measures of J = 



t(0)(/) 



By Lemma 5.7 and the following remark, the inequality 

<0(l){£i(J)/ 



t(o)(J) 



- £l 



holds with the minimal generalized i-characteristics {y'^{t;to,x)}xi£j of u'^ starting from J. 

In the i^-limit, by Ascoli-Arzelaa theorem the selected i^-characteristics {y'^{t;to,x)}xej converge to 
some generalized i-characteristics {y(t;to,x)}xe.i of u. By the upper continuity of Borel probability 
measures on compact sets one can derive 

-W(0)(J) =-lim<„„t(0)(j) 



< liniinf 0(l){/:i(J)/r + fii^ {A';'j) 
<0(l){£i(J)/r + /^^(X;)}. 



£1 



The same results hold for open sets by approximation from the interior. Then outer regularity yields the 
inequality for any Borel set. 

When to is one of the countably many points of 8 in Theorem |4. 2 [ one should take into account that 
t I— (wj)jump(O is continuous from the right, and thus ('yj:)cont- The thesis can thus be obtained by a 
diagonal argument for a sequence tk i to. □ 

Corollary 5.11. Consider a strictly hyperbolic system of conservation laws as and let the i-th 

characteristic field be genuinely non-linear, (2.1 : zGN). By Lemma 5.9 and the decay of positive waves 



in Theorem 10.3 of |Bre| . one has that there exists a positive constant C such that for every t > t > 0, 
every Borel subset BofM. and every solution u obtained as a limit of the front tracking approximation, 
the continuous part of the measures Vi satisfies 

\{v,Unt\{tm < CC\B)/T + C{^,','^'){[t 



W", 1^1 

(^)jump: (?/)coiit 
f 

A 

Xi{u) 

X[''\t,x), ~Kit,x) 
I, ff"^), r, f'-"^ 

CTl, . . . ,CTN, A, 

u 

u{t) 



T,t + T\ X Mj . 

Appendix A. Table of notations 

The Lebesgue measure of dimension respectively 1 and 2 

A positive constant which can be uniformly bounded in the limiting index 

Respectively positive part, negative part and variation of a measure t] 

The jump part and the continuous part of a measure rj 

Open subset of 

Flux function for a strictly hyperbolic system of conservation laws 
The Jacobian matrix A — Df 

The i-th eigenvalue oi A = D/, by strict hipcrbolicity Ai < • • • < An 
\i composed with u*^"^^ and a suitable mean at jump points (Page Island ( |4.3[ )) 
Left and right eigenvectors of A = D/, suitably normalized (Page pi 
Strengths of the waves and maps in the Riemann problem, see Page [i] and Q 
Mainly the variable in or the semigroup solution to 
The j/-front-tracking approximation, see recalls in Sect 
The (measure) derivatives of the BVioc(R^ x E) function u 
The restriction of the function u at time f > 0, see Page|4] 

Measures of the wave decomposition (2.5 1 of u^x^ along \ . . . , f^'' 



o ( O ) 
tionjij 



Conditional measures of v^^^ in the disintegration w.r.t. time, see e^. Pag ej? 
Parts of converging to the jump\continuous part of Vi (Section |4j Page 161 
The 'jump sets' of u and u'^ of Section |4j 



Interaction and interaction-cancellation measures, see (3.1) 
A w*-limit of ^il^ 

Balance and wave balance measures defined in Section 15.11 
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i^,out , (j ump\cont ) 



The interaction-cancellation-jump balance measures in (5.4 1, and a relative ?x;*-limit 



The fluxes introduced in the proof of Lemma 5.6 




The region bounded by selected (e.g. minimal) i-characteristics in (5.5) 
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